Computational investigations of the thermodynamic properties of size-selected water and Ar-water clusters: high-pressure transitions.

Classical parallel-tempering Monte Carlo simulations in the isothermal-isobaric ensemble were carried out for the (H2O)20 and Ar(H2O)20 clusters, over a wide range of temperatures (30-1000 K) and pressures (3 kPa-10 GPa) in order to study their thermodynamic properties and structural changes. The TIP4P/ice water model is employed for the water-water interactions, while both semiempirical and ab initio-based potentials are used to model the interaction between the rare-gas atoms and the water molecules. Temperature-pressure phase diagrams for these cluster systems were constructed by employing a two-dimensional multiple-histogram method. Structural changes were detected by analyzing the heat capacity landscape and the Pearson correlation coefficient profile for the interaction energy and volume. Those at high pressure correspond to solid-to-solid transitions and are found to be related to clathrate-like cages around the Ar atom. It is also shown that the formation and thermodynamic stability of such structures are determined by the intermolecular interaction between the rare-gas atoms and the host water molecules.

Weakly bound finite systems: (4He)N-Rb2(³Σu), clustering structures from a quantum Monte Carlo approach.

We report here ((4)He)(N)-Rb(2)((3)Σ(u)) complexes, 2 ≤ N ≤ 20, analysed through a quantum diffusion Monte Carlo stochastic approach. The calculations show that the spin stretched dimer molecule is bound outside the pure He sub-complex, due to the stronger He-He potential as compared with the He-Rb(2) interaction, while the rare gas atom moiety presents, in turn, a shell-like structure with ten He adatoms completing the first shell. Our results agree with previous findings on this and similarly weakly interacting systems.

Spin-polarized Rb2 interacting with bosonic He atoms: potential energy surface and quantum structures of small clusters.

A new full-dimension potential energy surface of the three-body He-Rb2(³Σ(u)(+)) complex and a quantum study of small (4He)(N)-Rb2(³Σ(u)(+)) clusters, 1 ≤ N ≤ 4, are presented. We have accurately fitted the ab initio points of the interaction to an analytical form and addressed the dopant's vibration, which is found to be negligible. A Variational approach and a Diffusion Monte Carlo technique have been applied to yield energy and geometric properties of the selected species. Our quantum structure calculations show a transition in the arrangements of the helium atoms from N = 2, where they tend to be separated across the diatomic bond, to N = 4, in which a closer packing of the rare gas particles is reached, guided by the dominance of the He-He potential over the weaker interaction of the latter adatoms with the doping dimer. The deepest well of the He-Rb2 interaction is placed at the T-shape configuration, a feature which causes the dopant to be located as parallel to the helium "minidroplet". Our results are shown to agree with previous findings on this and on similar systems.

Frequency domain description of Kohlrausch response through a pair of Havriliak-Negami-type functions: an analysis of functional proximity.

An approximation to the Fourier transform (FT) of the Kohlrausch function (stretched exponential) with shape parameter 0 < ß ≤ 1 is presented by using Havriliak-Negami-like functions. Mathematical expressions to fit their parameters α, γ, and τ, as functions of ß (0 < ß ≤ 1 and 1 < ß < 2) are given, which allows a quick identification in the frequency domain of the corresponding shape factor ß. Reconstruction via fast Fourier transform of frequency approximants to time domain are shown as good substitutes in short times though biased in long ones (increasing discrepancies as ß â†’ 1). The method is proposed as a template to commute time and frequency domains when analyzing complex data. Such a strategy facilitates intensive algorithmic search of parameters while adjusting the data of one or several Kohlrausch-Williams-Watts relaxations.

A new accurate and full dimensional potential energy surface of H5(+) based on a triatomics-in-molecules analytic functional form.

In this work a reliable full nine-dimensional potential energy surface for studying the dynamics of H(5)(+) is constructed, which is completely symmetric under any permutation of the nuclei. For this purpose, we develop a triatoms-in-molecules method as an extension of the more common diatoms-in-molecules one, which allows a very accurate description of the asymptotic regions by including correctly the charge-induced dipole and quadrupole interactions. Moreover, this treatment provides a semiquantitative description of all the topological features of the global potential compared with coupled cluster results. In particular, the hop of the proton between two H(2) fragments produces a double well in the potential. This resonant structure involving the five atoms produces a stabilization, lowering the barrier, and the triatoms-in-molecules yields to a barrier significantly higher than the ab initio results. Therefore, to improve the triatomics-in-molecules potential surface, two five-body terms are added, which are fitted to more than 110,000 coupled-cluster ab initio points. The global potential energy surface thus obtained in this work has an overall root mean square error of 0.079 kcal/mol for energies below 27 kcal/mol above the global well. The features of the potential are described and compared with previous available surfaces.

A path-integral Monte Carlo study of a small cluster: The Ar trimer.

The Ar(3) system has been studied between T=0 K and T=40 K by means of a path-integral Monte Carlo (PIMC) method. The behavior of the average energy in terms of the temperature has been explained by comparison with results obtained with the thermal averaged rovibrational spectra estimated via: (i) a quantum mechanical method based on distributed Gaussian functions for the interparticle distances and (ii) an analytical model which precisely accounts for the participation of the dissociative continua Ar(2)+Ar and Ar+Ar+Ar. Beyond T approximately 20 K, the system explores floppier configurations than the rigid equilateral geometry, as linear and Ar-Ar(2)-like arrangements, and fragmentates around T approximately 40 K. A careful investigation of the specific heat in terms of a confining radius in the PIMC calculation seems to discard a proper phase transition as in larger clusters, in apparent contradiction with previous reports of precise values for a liquid-gas transition. The onset of this noticeable change in the dynamics of the trimer occurs, however, at a remarkably low value of the temperature in comparison with Ar(n) systems formed with more Ar atoms. Quantum mechanical effects are found of relevance at T

An optimized full-configuration-interaction nuclear orbital approach to a "hard-core" interaction problem: application to (3He)(N)-Cl2(B) clusters (N < or = 4).

An efficient full-configuration-interaction nuclear orbital treatment has been recently developed as a benchmark quantum-chemistry-like method to calculate ground and excited "solvent" energies and wave functions in small doped DeltaE(est) clusters (N < or = 4) [M. P. de Lara-Castells, G. Delgado-Barrio, P. Villarreal, and A. O. Mitrushchenkov, J. Chem. Phys. 125, 221101 (2006)]. Additional methodological and computational details of the implementation, which uses an iterative Jacobi-Davidson diagonalization algorithm to properly address the inherent "hard-core" He-He interaction problem, are described here. The convergence of total energies, average pair He-He interaction energies, and relevant one- and two-body properties upon increasing the angular part of the one-particle basis set (expanded in spherical harmonics) has been analyzed, considering Cl(2) as the dopant and a semiempirical model (T-shaped) He-Cl(2)(B) potential. Converged results are used to analyze global energetic and structural aspects as well as the configuration makeup of the wave functions, associated with the ground and low-lying "solvent" excited states. Our study reveals that besides the fermionic nature of (3)He atoms, key roles in determining total binding energies and wave-function structures are played by the strong repulsive core of the He-He potential as well as its very weak attractive region, the most stable arrangement somehow departing from the one of N He atoms equally spaced on equatorial "ring" around the dopant. The present results for N = 4 fermions indicates the structural "pairing" of two (3)He atoms at opposite sides on a broad "belt" around the dopant, executing a sort of asymmetric umbrella motion. This pairing is a compromise between maximizing the (3)He-(3)He and the He-dopant attractions, and suppressing at the same time the "hard-core" repulsion. Although the He-He attractive interaction is rather weak, its contribution to the total energy is found to scale as a power of three and it thus increasingly affects the pair density distributions as the cluster grows in size.

A theoretical investigation on the spectrum of the Ar trimer for high rotational excitations.

A detailed study of the rovibrational spectrum of the Ar trimer is performed by means of an exact hyperspherical coordinate (HC) method and a variational approach based on distributed Gaussian functions (DGFs) to describe the interparticle distances. The good agreement observed between the energy levels obtained with both procedures for high values of the total angular momentum (J=15 and 20) reveals the quality of the DGF method to describe the rotation of the title system. Rotational constants for the lowest bound states, obtained as averages for each vibrational state, have been obtained and compared to previous results. A detailed analysis of density probability functions obtained by means of the HC approach for rovibrational states at J=0 and 20 shows close similitudes thus supporting the vibration-rotation separation adopted within the DGF scheme for the Ar(3) system.

Exact, Born-Oppenheimer, and quantum-chemistry-like calculations in helium clusters doped with light molecules: The He2N2(X) system.

Helium clusters doped with diatomic molecules, He(N)-BC, have been recently studied by means of a quantum-chemistry-like approach. The model treats He atoms as "electrons" and dopants as "nuclei" in standard electronic structure calculations. Due to the large mass difference between He atoms and electrons, and to the replacement of Coulomb interactions by intermolecular potentials, it is worth assessing up to what extent are the approximations involved in this model, i.e., decoupling of the BC rotation from the He-atom orbital angular momenta and Born-Oppenheimer separation of the BC stretch versus the He motions, accurate enough. These issues have been previously tackled elsewhere for the (4)He(2)-Br(2)(X) system, which contains a heavy dopant [Roncero et al., Int. J. Quantum Chem. 107, 2756 (2007)]. Here, we consider a similar cluster but with a much lighter dopant such as N(2)(X). Although the model does not provide the correct energy levels for the cluster, positions and intensities of the main detectable lines of the vibrotational Raman spectrum at low temperature are accurately reproduced.

A full-configuration interaction "nuclear orbital" method to study doped 3HeN clusters (N< or =4).

An efficient full configuration interaction (FCI) treatment, based on the Jacobi-Davidson algorithm, is developed in order to study small doped (3)He(N) clusters. The state of each He atom in a given cluster is described by a set of wave-functions which by extention of the quantum-chemistry notation are caller here "nuclear orbitals". The FCI treatment is applied to the calculation of binding energies and helium natural orbitals of (3)He(N)...Br(2)(X) complexes. In agreement with our previous calculations using a Hartree-Fock approach [Phys. Rev. Lett. 93, 053401 (2004)], in which the He-He interaction is modified at small distances to account for short-range correlation effects, the lowest-energy states of each multiplet are found to be very close in energy. The natural orbital analysis, in turn, indicates the adequacy of the "nuclear orbital" approach in these systems.

Vibrational quenching of CO2(010) by collisions with O(3P) at thermal energies: a quantum-mechanical study.

The CO(2)(010)-O((3)P) vibrational energy transfer (VET) efficiency is a key input to aeronomical models of the energy budget of the upper atmospheres of Earth, Venus, and Mars. This work addresses the physical mechanisms responsible for the high efficiency of the VET process at the thermal energies existing in the terrestrial upper atmosphere (150 K

The binding of 4He and 3He to a hydrogen molecule: a computational study for pH2 and oH2.

One bound state for the very weakly interacting complex of H(2) with one He atom is computed using a recently obtained, very accurate potential energy surface [A. I. Boothroyd and P. G. Martin, J. Chem. Phys. 119, 3187 (2003)] and employing a discrete variable representation expansion. The spatial features of this very diffused, "halo" state are analyzed and discussed for various isotopic variants of the title system.

A complete configurational study for the bound states of Ne trimers.

The structural properties and the energetics of the ground and the excited bound states of Ne(3) for zero total angular momentum are examined using different modelings for the two-body interactions. We employ a method consisting of a variational approach with a distributed Gaussian functions (DGF) basis set expansion. We discuss at length the advantages and possible limitations of such an approach, comparing it to other methods which have been applied in the literature to the same system. The DGF method turns out to be very accurate in giving us the bound states energetics and also provides in a natural way a convincing pictorial description of all the states, including those with dominant linear configurations. Additional bound states are found for the Ne(3) system with respect to those indicated in previous works and we suggest a "stabilization" procedure that can be used to assess the truly bound nature of a state. Some considerations on the relative reliability of the examined two-body interactions are also reported.

Bound-state energies in argon trimers via a variational expansion: the effects from many-body corrections.

In this paper we study the bound-state energies and geometries of Ar(3) for J=0, using the distributed Gaussian functions method that provides a configurational description of the different structures contributing to these states. Atom-atom potentials are employed and three-body long-range effects are also included in the computational treatment by adding to the sum of potentials the Axilrod-Teller triple-dipole correction for the whole rotationless energy spectrum. An estimate of the total number of bound states for the Ar trimer is given. With respect to previous calculations, limited to the lower-lying states, our results show slightly larger nonadditive effects and are further able to predict the full range of the bound spectrum. Changes on the geometries of a large part of the vibrationally excited states of Ar(3) when the Axilrod-Teller term is included in the molecular potential are found by the present study.

Role of boson-fermion statistics on the Raman spectra of Br2(X) in helium clusters.

The role played by the bosonic or fermionic character of He atoms surrounding a Br2(X) molecule is analyzed through vibrotational Raman spectra simulations. Quantum chemistry-type calculations reveal the spin multiplicity to be chiefly responsible for the drastic difference observed by Grebenev et al. [Science 279, 2083 (1998)]] in the rotational structure of molecules embedded in helium droplets.

Raman spectra of (He)N-Br2(X) clusters: the role of boson/fermion statistics in a quantum solvent.

The aim of this paper is to elucidate the role played by the bosonic/fermionic character of N He atoms solvating a Br2(X) molecule. To this end, an adiabatic model in the molecular stretching coordinate is assumed and the ground energy levels of the complexes are searched by means of Hartree (or Hartree-Fock) Quantum Chemistry calculations for 4He (or 3He) solvent atoms. Simulations of vib-rotational Raman spectra point at the spin multiplicity as the main feature responsible for the drastic difference in the rotational structures of molecules embedded in boson or fermion helium drops as already observed by the experiments of Grebenev et al. [S. Grebenev, J. P. Toennies, and A. F. Vilesov, Science 279 (1998) 2083].

The open-shell interaction of He with the B 3Piu(0+) state of Br2: an ab initio study and its comparison with a diatomics-in-molecule perturbation model.

The interaction of He with Br2 in electronically excited B 3Piu state is investigated using spin-unrestricted single and double coupled-cluster approach with noniterative perturbative treatment of triple excitations. Internal electrons of the Br atom are described by effective core pseudopotentials. The validity of this approach is analyzed by comparing the lowest 2Sigma+ and 2Pi electronic states of the HeBr molecule with those obtained in all electron calculations [J. Chem. Phys. 115, 10438 (2001)]. In this context, we examine the performance of different basis sets and saturation with bond functions. The comparison of theoretical blue-shifts with the experiment provides confidence about the present ab initio calculations. In addition, He-Br results of ab initio calculations at the same level are used to obtain approximate He-Br2 (3Piu) interactions in the framework of the diatomics-in-molecule first order perturbation theory (IDIM-PT1) [J. Chem. Phys. 104, 9913 (1996)]. Overall, the IDIM-PT1 model results show a good agreement with the ab initio ones, being the main difference the sensitivity to the elongation of the Br-Br bond.

Matrix elements for anharmonic potentials: Application to I2 morse oscillator.

Many authors have contributed expressions for obtaining analytical matrix elements for Morse oscillators. In this work, we discuss the advantages of using these expressions. At the same time, we propose a full numerical method to calculate these matrix elements and we compare, for the I2 system, the different results given by Gallas, Vasan, and Cross, and the variational method.